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ABSTRACT 


iii . 


The  bound  state  of  two  massive  boson 
fields  coupled  together  by  a  real  scalar  field  is 
studied.  Using  the  MN-quanturn  approximation",  an 
integral  equation  for  the  "wave  function"  is  obtained. 
The  eigenvalues  of  the  resulting  equation  are  com¬ 
pared  with  the  eigenvalues  of  the  Schrodinger,  Klein- 
Gordon,  Bethe-Salpeter  and  Son-Sucher  equations. 

'When  the  binding  energy  is  small  compared  with  the 
masses  of  the  fields,  the  coupling  constant  obtained 
from  this  equation  is  near  to  those  obtained  from  the 
Klein-Gordon  and  Schrodinger  equations,  but  is  smaller 
than  the  coupling  constant  obtained  from  the  Bethe- 
Salpeter  and  Son-Sucher  equations. 
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NOTATIONS 


Throughout  this  work  we  shall 
=  {1,-1, -1,-1}  .  We  shall  also  use 

notations ; 

x,k,k*  -  4-vector 

4  4  4 

dx  ,  dk  ,  dk'  ,  to  be  d  x  ,  d  k  ,  d  k' 
p  j  P1  j  q  ^  r  ,  r'  -  3-vector,  with 
rx  to  be  their  magnitude  respectively. 


-  E  -  The  binding  energy 
2 

G  /4t t  -  The  coupling  constant 

2  2 
i  +  p  +  q 

2pq 

0,  2,  2x1/2/  2,  2x1/2  ,  n  0  2 

2(p  +m2)  (q  +m2)  +  1  -  2m2 

— 


B  = 
A  = 

t  = 

t  = 


use  the  metric 
the  following 

respectively 

P  ,  Pf  .  q  *  r  , 

Also 


-  ■ 


CHAPTER  I. 


INTRODUCTION 


The  binding  energy  of  a  system  has  always 
been  a  quantity  of  great  interest  to  the  theorists  in 
physics,  since  it  could  usually  be  measured  experi¬ 
mentally.  Thus,  any  theory  attempting  to  explain  the 
bound  states  could  be  immediately  tested  by  experiment 

With  the  introduction  of  quantum  mechanics, 
the  binding  energies  of  electrons  in  different  systems 
such  as  hydrogen  and  helium  could  be  explained  reason¬ 
ably  well  using  the  Schrodinger  equation  with  the  cor¬ 
rect  form  of  the  potential''".  Following  this  success, 
people  have  tried  to  extend  this  method  to  systems 
involving  nuclear  interactions.  Hence,  solutions  of 
the  Schrodinger  equation  using  simple  potentials  such 
as  the  square  well,  exponential  and  gaussian  potential 

were  used  to  explain  nuclear  bound  state  systems  such 

2 

as  the  deutron  .  This  method  had  two  important  draw¬ 
backs.  First,  to  justify  the  use  of  the  Schrodinger 
equation  one  had  to  assume  that  the  bound  particle  had 
a  nonrelativistic  velocity;  and  second,  one  had  no 
theoretical  grounds  to  assume  that  any  of  these  potent 
ials  actually  represented  the  nuclear  forces . 


. 


In  1935*  Yukawa-3 ,  in  trying  to  formulate  a 
consistent  relativistic  theory,  had  to  assume  the 
existence  of  a  particle,  now  called  a  meson,  having  a 
rest  mass  between  the  electron  and  nucleon  masses. 

He  proposed  a  theory  of  nuclear  forces  which  involves 
these  mesons,  and  is  referred  to  as  the  "meson  theory 
of  nuclear  forces".  The  detection  of  such  a  particle 
in  the  cosmic  radiation  and  the  subsequent  discovery 
that  there  exist  several  different  mesons  with  different 
charges  and  rest  masses,  stimulated  many  attempts  to 
fit  these  mesons  into  a  consistent  scheme  of  nuclear 
forces.  These  attempts,  however,  have  not  yet  succeed¬ 
ed  in  reproducing  quantitatively  the  known  properties 
of  nuclear  forces.  Thus,  while  there  is  still  no  gen¬ 
eral  theory  for  the  interaction  of  elementary  particles, 
people  have  recently  proposed  different  relativistic 
equations,  based  on  different  assumptions  or  methods, 

which  still  have  no  solid  theoretical  foundation.  Two 

4 

such  equations  are  the  Bethe-Salpeter  equation  and 

5 

a  bound  state  equation  proposed  by  Son  and  Sucher  . 


In  a  paper  by  Greenberg 


6 


a  new  method  has 


been  developed  to  obtain  approximate  solutions  for  the 
equations  of  motion  of  quantized  fields.  This  tech- 


9 1 

' 


3 


nique,  the  "H'-quantum  approximation"  is  a  covariant 
method  of  field  quantization  which  is  suitable  for 
strong  interaction  physics,  where  the  results  of 
perturbative  calculations  are  not  reliable. 

This  work  is  an  application . of  the  "N- 
quantum  approximation".  We  propose  to  determine  the 
binding  energy  of  two  massive  spinless  fields  and 
^2  coupled  together  by  a  real  scalar  field  $  of 
unit  mass.  A  similar  problem  for  two  spinor  fields 
coupled  together  by  a  pseudoscalar  meson  field  is 

7 

discussed  in  a  paper  by  Greenberg  and  Genolio1 .  We 
shall  derive  the  equation  for  the  "wave  function" 
of  the  bound  system  and  refer  the  reader  to  refer¬ 
ences  (6)  and  (7)  for  a  careful  discussion  about  the 
assumptions,  the  nature,  and  the  range  of  validity  of 
the  approximation. 

For  the  sake  of  completeness,  we  outline 
here  the  essential  features  of  the  method.  We  assume 
that  we ' are  dealing  with  a  field  theory  where  there  are 
irreducible  sets  of  field  operators.  The  irreducible 
sets  of  fields  are  the  Heisenberg  fields,  which  appear 
in  the  Lagrangian  and  Hamiltonian,  but  which  have  no 


. 


4. 


simple  universal  relation  to  the  particles*  states; 
and  the  In  (out)  fields,  which  do  not  appear  in  the 
Lagrangian  or  the  Hamiltonian,  but  are  the  free 
fields  describing  the  asymptotic  in  (out)  states  of 
freely  moving  particles.  The  basic  idea  of  the  "N- 
quantum  approximation"  is  to  exploit  the  expansion  of 
the  Heisenberg  fields  in  terms  of  normal -ordered 

o 

products  of  in-fields  as  was  first  given  by  Haag  . 

For  an  exact  Heisenberg  field  of  any  non-trivial  loc¬ 
al  field  theory,  this  expansion  never  terminates;  none¬ 
theless,  the  expansion  can  terminate  for  an  approximation 
to  the  exact  Heisenberg  field.  Thus  by  substituting 
the  truncated  expansions  for  the  Heisenberg  fields  into 
their  equations  of  motion  and  using  the  fact  that 
various  normal  ordered  products  of  in-fields  are  in¬ 
dependent,  one  obtains  equations  for  the  expansion 
coefficients.  In  general  these  equations  are  coupled 
sets  of  nonlinear  singular  integral  equations.  For 
the  lowest  order  approximation,  however,  one  obtains  a 
linear  equation. 

The  advantages  of  this  approach  are  that  the 
resultant  equations  are  separable  and  relatively  simple; 
that  is  to  say,  the  time  component  can  be  separated  out. 


. 


5. 


This  is  not  the  case  for  the  Bethe-Saipeter  equation. 
This  method  is  also  non -perturbative *  and  therefore  it 
does  not  depend  on  an  expansion  in  terms  of  the  large 
coupling  constant.  Also  higher  order  correction  can 
be  added  in  a  straightforward  manner  as  indicated  in 
references  (6)  and  (7) • 

In  this  work  we  determined  the  binding 
energies  predicted  by  the  lowest  order  approximation* 
and  compared  them  with  those  obtained  from  the  solutions 
of  Schrodinger*  Klein-Gordon*  Bethe-Saipeter  and  Son- 
Sucher  equations. 


Chapter  II.  SCHRODINGER  AND  KLEIN -GORDON  EQUATIONS 

IN  MOMENTUM  SPACE 


We  transform  Schrodinger  and  Klein-Gordon 
equations  into  momentum  space  and  obtain  two  linear 
integral  equations  which  will  be  solved  numerically  In 
chapter  IV. 

II. 1  Schrodinger  Equation 


Schrodinger  equation  for  two  particles  of 
mass  and  in  the  relative  coordinates  is 

-  ^  V2f(r)  +  v(r)?(?)  =  E?(r)  (Il-l) 

2 

with  ft  =  1  and  M  =  m-^m^/ (m-^+rn^)  . 

We  define  the  Fourier  Transform  of  u(p)  and  v(p’)  by 

_3 

f(r)  =  (2 it)2  /  e“ip,ru(p)d^p  =  F[u]  , 

_3 

v (r)  =  (2tt)2  /  e_1P-r  v(p)d3p  =  F[v]  . 


and 


Substituting  these  definitions  into  equation  (il-l) 
we  get 


r  2 

P_ 

2M 


u(p) 


+  F 


(2tt)2  /  vfp’-p1  )u(p 1 )  d^p  * 


=  EF[u(p)]  * 


where  we  have  used  the  convolution  theorem 


3 


F[A]F[B]  =  F 


(2t r)  /  A(p-p  '  )B(pl  )d'3p 


t  'l  t 


Taking  the  inverse  Fourier  Transform  and  letting 
E  =  -B  where  B  >  0  *  we  obtain 


3 


[§M  +  bMp)  =  -(27t)^  /  v(p-p*  )u(pT  )d'3pt  * 


2M 


t  \  a^-r,  t 


(II-2) 


where  B  is  the  binding  energy  of  the  system*  and 
u(p)  is  the  bound  state  wave  function.  Since  v(p-pl ) 
is  invariant  under  rotation  we  may  write*  with 
x  =  p.p'/pp1  , 


00 


V 


(p-pr)  =  2  (2/e+l)v^(p*pt  )P  p(x)  * 


i=0 


i 


(H-3) 


where  ( x )  are  the  Legendre  polynomials.  Using  the 


identity 


. 


8. 


1 


f  P£t(x)Pj,(x)dx  =  2I+I 


-1 


(II -4) 


we  obtain 


v,(p,pT)  =  i  /  v(p-pl )P.  (x)dx 


4 


-1 


£ 


(H-5) 


We  now  set 


u(p)  =  s  u^(p)Y^m(p) 


(n-6) 


where  Y.(p)  are  the  usual  spherical  harmonics. 
Substituting  equations  ( II —3 )  and  (II-6)  into  equation 
(II-2),  and  using 


/  ee,(p.$')YeJp')<ip'  =  W'pU 


£m 


?JW1  £mKy/u££ 


(II-7) 


we  get  the  reduced  equation 


2  —  o° 

[fjj  +  B]u^(p)  =  -(f)2  /  v^(p,p')u^(p')P'  dp' 


(XI -8) 


For  spherically  symmetric  potentials. 


v ( p )  =  (27 r)?  /  elr’*p  v(?)d^r 


3. 


\ 


Q 


-i  oo  1  . 

=  (2ir)  f  r  v(r)dr  f  e 1Pr^  by 
0  -1 

1 

2  2  °°  -1 

~  (— )  /  P  sin(pr) v(r)rdr  . 


0 


Therefore,, 


1 


/-+•  t,\  /2n2  r°°  sin  |p-pl|  r  ,  \  , 

v(p-pT)  =  (-)  / - — —  v(r)  rdr 

0  lp-p  'I 


(II -9) 


Using  the  identity 


00 


Sj-..n..[p.  P'l  r  =  TT  2  (pp> ) 2 (m-t-l/2 )  /0(pr)  J  ,Q(p'r) 


IP-P 


m=0  Wl/2'*“-'VV2 


PmW  * 


where  ^m+l/2^z^  ^esse^-  function,,  and  equations 

(II-9)  and  (II.-4),  we  can  write  equation  (II-5)  as 


vi(p,pt)  =  (f)^  /  (PP T )  f+l/2  (p  ' r) v  (r)  rc^r 


1  1 
—  00  _ — 


(11-10) 


Equation  (II-8),  with  v^(p,pT)  given  by  equation 
(II  -10) j  is  our  Schrodinger  equation  in  momentum  space 
and  we  shall  solve  it  later  on  for  different  potentials 
v(r)  . 


10. 


II. 2  Klein-Gordon  Equation 

We  write  the  Klein-Gordon  equation  for  a 
particle  of  mass  m  ,  in  a  central  field  of  force 
v(r)  j  as 

[A  +  m^]¥(r,  t)  =  -2mv(r)3r('rJt)  ^  (il-ll) 

where  the  factor  2m  is  necessary  to  give  the  correct 
non-relativistic  limit  of  equation  (II-2) .  Let 
5f(r,t)  =  f( r)e  , 

where  E  =  m-B  is  the  bound  state  energy.  Then 
equation  (il-ll)  becomes 

V^V(r)  +  (E2-m2)^(r)  -  2mv(r)^(r)  =  0  , 


or 

-  ^  V>(r)  +v(r)^(r)  =  -  f  (?)  .  (H-12) 


Note  that  equation  (11-12)  is  identical  with  equation 
(Il-l)  provided  we  replace  in  equation  (il-l)  -E  =  B 
by  B(2m-B)/2m  and  2M  by  2m  .  Hence,  the  analogue 


. 


11. 


of  equation  (II-8)  is  now 


r ,  B(2m-B)n  f  \ 

[35  +  2i-^ufP)  = 


1 


00 


-(§)2  /  v((p,p’)u((p')p'2dp' 

(H-13) 


'T  J0  !~i 


where  v^P-jP1)  is  still  given  by  equation  (11-10) 


II. 3  The  Yukawa  Potential 


Consider  the  Yukawa  potential 


v  (r)  = 


p“!“Lr 

-A—  , 

\i  r 


where  A  =  is  the  coupling  constant  and  we  choose 

p  to  be  unity.  Equation  (II-IO)  now  becomes 


1 


CO 


1 


vx(p,PT)  =  (PP  1 )  2J  Jl+l/2 


0 


i+l/2  (Pr)  tJi+l/2 (p  T  r)  e_rdr 


1 


p-TT-  N  ^  ^  n  fl+P  +P*  > 

'  1  ppl  uP  2pp f  ' 


where  Q^(z)  are  the  Legendre  functions  of  the 
second  kind.  Thus  equations  (II-8)  and  (11-13)  now 


become 


_ 


12. 


2  00 

[fpj  +  B]Ug  (p)  =  |  /  Q^(t)u^(q)q  |S.  ,  (II-14) 

and 

tfs +  ^sr“]u2p)  -  ?  /0Qi(t)ue(q)q  f1  (II"15) 

respectively,  where 

-,,22 

4.  _  1+P  +q 

r  ~  2pq  ~  * 


13. 


Chapter  III.  NEW  RELATIVISTIC  EQUATION 

FOR  BOUND  STATE 

In  this  chapter  we  shall  derive  a  relativ¬ 
istic  equation  using  the  lowest -order  MN -quantum 
approximation'’  .  This  equation  reduces  to  the  Schrodin- 
ger  and  Klein-Gordon  equations  in  the  proper  limits. 


III.l  Equations  of  Motion 

Let  the  Hamiltonian  for  two  scalar  fields 
and  masses  and  respectively,, 

which  interact  with  the  exchange  of  a  scalar  meson  0 
of  unit  mass,  be: 

H  =  /  H  dx 


2 

=  f  Z  [IT,  (x )  TT .  (x)  + 
i=l 


(VTi)(VTi)  +  m^Ti(x)Ti(x) ]dx 


+  4  /[ P2(x)  +  $2(x)  +  (V$)2]dr  +  G  f  Z  m.  )dT, 

d  i=1  ±  x  ± 

(III-1) 


14. 


where  H  is  the  Hamiltonian  density.  From  Hamilton’s 
canonical  equations,  it  follows  that 


1 

3t 


dH 

W. 


3 


_  Z 


c3x  .  (W . 


=  7 r. 


and 


c>t 


M 

c)f. 


3  . 

V  O 

,_1^. 

J-1  J 


dH 


=  m25h  +  2Gmi$T1  -  V2^ 


Therefore ^ 


St2 


m?f.  +  2Gm.  -  V2F. 

li  l  1  i 


or. 


15 


(HI -2) 


Similarly,  the  equation  of  motion  of  the  meson  field 
turns  out  to  be 


2 

[  A  +  1  ]  $(x)  =  -  G  2m.[!.,?.] 

i=l  111 


(Ill -3) 


Equations  (III-2)  and  (III-3)  are  the  equations  of 
motion  of  the  fields.  Note  however,  that  equation 
(III-2)  is  not  symmetric  in  $  and  55  ,  and  since, 

from  equation  (ill -3)*  $  will  depend  in  general 

upon  35  and  'L  ,  we  would  like  to  symmetrize 
equation  (ill -2)  to  prevent  any  ambiguity.  Since 
$  and  35  commute,  we  can  rewrite  equation  (III-2) 


as 


[  A  +  m?  ]  3/i(x)  =  -  Gnu  {$,35}  . 


(III-4) 


Equations  (III-3)  and  (ill -4)  can  be 


■  . 


1 6. 


transformed,  into  momentum  space  by  defining 
<t>(x)  =  /  O'  (k)e~lkx  dk  , 

^(x)  =  /  3'!(k)e'lkx  dk 

and 

f±(x)  =  /  H (k)e~ikx  dk  . 

Equation  (III-4)  then  becomes: 

/5'1(k)  (mhk2)e“lkxdk  =  -  am±/  [O^ ) ,  (kg)  ]e"^kl+k2)x 

dk1dk2 


where  we  have  dropped  the  primes  on  $  and 

Multiplying  both  sides  of  the  above  equation  by 
i  k1  x 

e  dx  ,  integrating  over  x^  and  using  the  identity 
f  )x  c[x  =  (2ir) ^6  (k-k* )  ; 

we  get 

(m2-k2)f1(k)  =  -  Gm1/B(k-k1-k2)  ($(k1),l'1(k2))dk1dk2 


(III-5) 


' 


17. 


Similarly  equation  (ill -3)  becomes 


(l-k2)<t>(k)  = 


-  G  2  m1  ;  Sfk-^-kg)  (i'i(k1),f1(k2)  Jd^dkg  .  (III-6) 

Equations  (ill -5)  and  (ill -6)  are  the  equations  of 
motion  of  the  Heisenberg  fields  in  momentum  space. 


Ill .2  The  "N-Quantum  Approximation 


Now  we  want  to  write  down  the  truncated 
expansion  for  the  Heisenberg  fields  in  terms  of  the 
in-fields.  Since  3^  ,  ,  anc^  ^  rePresent  particles 

which  are  stable;  the  set  of  in-fields  will  include 
,  and  <£>^n  .  Each  stable  bound  state 

will  also  lead  to  an  in-field,  and  since  we  are 
dealing  with  the  lowest -order  "N-quantum  approximation", 
which  neglects  any  radiative  transition,  all  our 
bound  states  are  stable.  Hence,  the  set  of  in-fields 
will  include  ,  the  bound  state  in-fields,  where 

n  is  the  principle  quantum  number  and  l  is  the 
angular  momentum. 


. 


18. 


To  obtain  the  lowest -order-  "N-quantum  approximation" 

X  y 

we  use  following  Greenberg rs  work  ' ' ,  the  following 
truncated  expansions; 


4>(k)  =  4>ln(k)61(k) 


+  /d^dkgdfk-^-kg)  Z  g1(k1Jk2):^n(k  )6  (n)?T(k2)6  (k  )  : 

i_l  iUi  -1-  ^  lui 


and 


(III-7) 


?1(k)  =  <t'P(k)6m  (k)  +  /dkjdkgSCk-k-^kp  2  fn/k:,k2) 

1  n,  i 

:4n(kf  V<k l)DnIV5d  /k2)  =  .  (III-8) 

c.  n  z 

where  f  .(k-,,^)  and  g.(kn,kn)  are  the  vertex  functions 
niv  1*  2'  °iv  1  27 

which  are  related,  in  the  weak  coupling  limit,  to  the 
matrix  elements 


T-j_  PnXi>  and  respectively. 


in 


.in,  ,  =in, 


2  ,_2x  ^in 


Also  6m(k)  -  6(m  -k  )  ,  are  the  bound  states1 

in-fields,  and  d^  =  £  are  masses  the 

bound  states  with  the  binding  energies  B  ^  .  There 

are  similar  expansions  for  and  5h  but  we  shall  not 

need  them  here.  In  what  follows  we  shall  often  abbreviate 

?Jn(kn)5  (k,  )Dll?(k0)6  ,  (k„)  to 
2  v  1 ’  m^  1  n^v  2'  dfi^  2' 


an  expression  such  as 


19. 

^2 )Dn£ when  there  is  no  danger  of  an  ambiguity. 
Substitution  of  equations  (III-7)  and  (III-8)  into 
equation  (III-5)  yields: 

/dk1dk2(m®-k2)6(k-k1-k2)  Z  fni(k1,k2):?2(l)Dni(2): 

=  -  Gm1/dk1dk26  (k-k1-k2)  {  O(l)  +  /dk^dk^b  (k-^k  -k^) 

2 

*  .21si^k3"k4^  :yi(3)^i(4)  :  ] ,  [^(2)  +  /dk^dk^b  (k2-k5~k^) 

fn£(k5^k6)  :f2(5)Dni(6)  :]  ]  *  (HI -9) 

Since  normal  products  of  in-fields  are  independent  of 
each  other,  we  shall  try  to  equate  the  coefficient  of 
:f2 (l)Dn^ (2)  :  on  the  left  hand  side  of  equation  (III-9), 
with  those  coefficients  of  the  normal  products,  on  the 
right  hand  side  of  equation  (111-9),  which  are  of  the 
same  order.  Since  we  are  interested  in  the  lowest-order 
approximation  we  shall  drop  any  higher  terms  of  the 
normal  products.  Hence,  the  only  term  to  survive,  from 
the  right  hand  side  of  equation  (III-9)  is 
C  :'i'1(3)^1  (4)  +  i'2(3)?2(4)  :f2(5)Dn^(6)  :]  .  The  term 

( :T-,  (3) q  (4 )  : ,  :?2 (5)D  » (6)  : )  however,  cannot  be  further 


contracted  to  lower  order  and  thus,  in  the  lowest-order 


- 
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approximation,,  it  will  be  ignored.  Therefore,,  we 
are  left  with 


t:S'2(3)?2(4):J:?2(5)D  ni(6):} 


=  :?2(3)?2(4):?2(5)D ^(6):  +:?2(5)D^  (6)  :'F2(3)f2(4)  : 

=  V2(3)?2(5):?2(4)Dae(6):  +  ^(5)?2(3)  :52(4)Dni(6)  :  . 

(Ill -10) 


To  find  the  explicit  form  of  the  contraction  ¥0(3)^2^) 

c - 1 

and  we  use  the  commutation  relation; 


[fin(k1)6m(ki),fln(k 


2>Vk; 


)3  =  (2ir) 


-3 


s(ki)Sm(ki)5(ki+k2) 


(HI -11) 


+1  for  k<^  >  0 
where  e(k)  =  C  n 

-1  for  ku  <  0 


Talcing  the  vacuum  amplitude  we  get 

<fn(ki)8m(ki)fn(k2)6m(k2)>o=  (2Tr)-3S(k1)8m(k1)6(k1+k2) 


(III-12) 
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1  for  k°  >  0 

where  0(k)  =  {  n  .  Then, 

0  for  ku  <  0 

<fln(k2)5m(k2)'i'ln(ki)6m(ki)>£>=  (27r)-3S(k2)6m(k2)6(k1+k2) 

(III -13) 

and,  as  can  be  verified  directly,  equations  (III-12) 
and  (III-13)  satisfy  the  commutation  relation  (ill-ll) . 
Therefore,  using  equations  (III-12)  and  (III-13)  we  can 
write  equation  (III-IO)  as 

C  :9'2(3)$2(4)  ■V2(5)l>nJl(6)  :) 

=  (2ir)'36  (3)5(k3+k5)[0(k3)+0(-k3)]  :f2(4)Dny6) 

=  (2tt)’36  (3)8(k3+k5)  :?2(4)Dni(6)  :  .  (III-14) 

Equation  (ill -9)  now  becomes 

/d^dkgOt^-kb^k-k-L-ky  2  (k^kg)  :?2(l)Dn<e  (2)  : 

n,  i 

Gm, 

=  -  — — ^/d^  .  .  .dk^6  (k-^-kg)  6  (k1-k^-k^)  6  (k2-k^-k 

:k^4)Dni(6)  : 

n ,  Jo 
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where  we  used  the  fact  that  we  can  interchange  the 
order  of  integration  and  thus,  relable  the  dummy  variables 
k^-k^  as  k^  ,  kg+k^  as  k2  ,  and  k^  as  -k^ 

Since  equation  (III-15)  is  true  for  all  k  ,  n  ,  and 
i  we  then  have; 

p  p  Gmi 

U1-(k1+k2)  ]fnJ^(k13k2)  =  "  '^y3^dk3fni  ^k3"k2^S2^~k3"kl^  6m2^k3^ 

(III -l6) 

2  2  2  2 

with  the  constrains;  k1  =  m2  and  k2  =  dn^  .  To 

calculate  g2(-k^,k^)  explicitly,  we  substitute  equations 

(III -7)  and  (III-8)  into  equation  (III-6),  and  retain 

only  the  first  term  on  the  right  hand  side  of  the  expansions 

for  and  'L  ,  which  amount  to  just  the  Born  approximation. 

Thus, 


'  . 
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p 


/dk1dk2(l-k  ) 6  (k-k-j  -kQ)  2  g,  (kn  , kQ)  :  i,  (l)  i.  (2)  : 


i  2 7  i!Z1  i  '1*  2 


G  2  m^^/d^dkgB  (k-k-^k^  :  { 5k  (l)  ,  5^  (2) }  : 


(HI-17) 


Let  us  write  i^n  and  ifn  in  terms  of  the  creation  and 
annhilation  operators ,  it  and  5k  respectively,  as 


^n(ki)  =  N[yt(l)  +  5k  ( 1 )  ]  , 


■i  v  1 


i 


and 


^n(k2)  =  N[i'1(2)  +  f(2)] 

where  N  is  the  normalization  constant.  Taking  the  vacuum 
state  expectation  value  of  both  sides  of  equation  (111-17), 
and  using  the  properties  of  the  creation  and  annhilation 
operators,  and  the  orthogonality  property  of  states j  we 
get 


[1  -  (k1+k2)2]gi(k1,kp)  =  -  2Gmi 


or 


si  (k^kp) 


2Gm. 

_ i 

1  -  (kx+k2)2 


(III -18) 
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Thus,  equation  (III-16)  becomes 


[mi  '  = 


(HI-19) 


We  require  that  equation  (III-19)  should  have  the 
correct  non-relativistic  limit.  This  can  be  achieved 


by  choosing  k2  >  0  and  k^  <  0  .  Hence,  we  rewrite 

equation  (III-19)  as 


fni  k2^ 


-  (k1+k2)2]fn^(k1,k2)  = 


1  -  (k-,-k)2 


(III-20) 


To  see  that  the  factor  2  on  the  right  hand  side  of 
equation  (III-20)  is  the  correct  one,  we  note  that  if 
we  have  used  the  unsymmetrized  form  of  equation  (III-4), 
namely,  equation  (III-2);  we  would  have  obtained  equation 
(III-20)  directly. 


The  conditions 


25. 


suggest  that  we  choose  our  coordinates,  as 


1 


=  f-  (p2+»|)2i  p) 


k2  =  C  dni  ;  ^  ] 


i  r  /  2  2x2  — ►  ■> 
k  =  l  -(q  +mQ )  ;  q  J 


Then,  equation  (ill -20)  reduces  to 


[m. 


-  m 


1 

-dh  +  2dni(p2+tr|)2]f„s(p) 


2  ni 


ni 


4G^m, mQ 

— — /dk6m  (k)0(-k)- 
(2tt)j  m2  1 


Wk) 


m2  +  2k^k 


-  m. 


With 


1 


1 


2k1k 


2(k°k°  - 


2  ( p2+m2 )  2  ( q2+m2 )  2  -  2]?1.q 


and 


6  (k) 

m2  v  ' 


K  /  2  .  2x  _  ,  2  .  02^2n 

6(mQ-k  )  =  5(mQ-k  +k  ) 


1  1 

e  f  r  /  2  2x2  ,  0]  r  /  2  2x2  ,  0n 

6  [  [  ( q  +mQ )  -  k  J  L  ( q  +mQ )  +  k  J  J 


4(q2+m2)2{6[k°  -  (q2+m2)2]  +  5[k°  +  (q2+m2)2]  ) 


. 
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we  have 


td ni  +  ‘  2d^(pP+m2)2]f^(p) 


1 


ni v  ^  2 ■ 

1 


'ni 


1 


2 

-  2Q  m-,  mQ  q  q  o — o  2  n 
:  - i2  /d3q(q2+m2)2  /  dk° 

(2tt)3  2 


6[k°  +  (q2+m2)2]  fn^(q) 

I  I 

1  -  2m2  +  2  ( p2-Hn2 ) 2  ( q2+m2 ) 2  -  2p.q 
1 


2 

2G  mn m0  0 

1  2  fa3 q 


(q2+m2)^  fn^(q) 


(2tt) 


3 


1  1 

„  2  0  ,  2  2\2/2  2\2  0-»  -> 
1  -  2m2  +  2 (p  +m2)  (q  +m2)  -  2p.q 


(III -21) 


Let 


1 


1 


D  =  1  -  2m?  +  2  (p2+rru) 2  (q2+m?) 2  -  2p.q 


=  2pq[ 


*2'  vh  '-2' 

1  1 
0  2  ,  o  /  2 ,  2x2,  2,  2x2 

1  -  2m2  +  2 ( p  +m0 )  ( q  +m2 ) 

2pq  T 


P.q 

pq 


=  2pq(t-s) 


where 

1 

t  =  — 


1  1 
0  2  l0,  2,  2x2,  2,  2x2 

2m2  +2(p  +m2)  (q  +m2) 
2pq 


and 
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s  S  1  . 

pq  - 


Using  the  identity 


00 


t  -  S 


=  2  (2i+l)Pi(s)Q^(t) 


£=0 


<  1 


we  get 


00 


b  =  2k^0(2i+1)Vs)Ie^ 


Hence,  equation  (III-21)  reduces  to 


1 

[dni  +  m2  -  ml  -  2dnUp2+m2)2]fnpP) 


n£ 


2  ' 


G  mi mQ  °°  o  Q;(t)  f  p  (q) 

2  P  (s)(2J+l)  / d3q  -i - 


(2ir)3  ^=0 


/  2,  2N  2 

pq(q  +m2) 


Let 


fnPp)  =  UHp3W^ 


using  equations  (II-4)  and  (II-7)  we  finally  obtain 


1 

[dL  +  m2  -  ml  -  2dnLp2+m2)23fnLp) 


°2T2t  qV^)fn,e(q) 

- - — j  dq - 


2t t  0 


1 

(  2,  2^2 
p ( q  +m2 ) 


(ill -22) 


. 


. 
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where  we  have  dropped  the  prime  on  f  .  Equation 
(ill -22)  is  our  relativistic  equation  and  we  shall 
solve  it  numerically  later  on. 


III. 3  Reduction  to  Schrodinger  Equation 

Consider  equation  (III-22)  in  the 

non-relativistic  limit  when  the  masses  of  the  particles 

are  much  greater  than  their  momenta  and  binding 

energies.  Then, 

1  1 

0  /  2  2x2/  2  2x2  n  0  2 

2 (p  +m2)  (q  +m2)  +  1  -  2m2 

-  =  2pq 


2pq 


[  2m2(l  +  -Zg-)  (1  +  -\)  +  1  -  2m2  ] 


1 


2m 

P 


2  2m2 
2  2 


2 nri  ^  Pm2  (  q  +  ~^2  +  +  1  -  2m2  ] 

dpq  d  2m2  2m2 


2  2 

_  p  +  q  +  1 
2pq 


t 


1 

/  2  2x2 

(q  +  mp) 


1 

m. 


(1  - 


2rrL 


1_ 

m. 


and 
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4  +  m2  -  ml  - 


2.2  2 


1 

2.  2x2 


=  (m1+m2-Bni)  +  m2  -  m1  -  2 (rm^-B^)  (p  +m*) 


2m2  +  2m1m2  2Bn^mi  2Bnim2  +  Bn &  ~  (mi+m2  Bni^2m2 


2 

X_) 

iru' 


2  m  +mp 

2m,  (  §  — - -  +  B  g  )  = 

1 x  2  m1m2  nH  ' 


2m1^2M  +  BnP 


Thus,  in  the  non-relativistic  limit,  equation  (ill -22) 
reduces  to 


n2  oo 

G  If  j  _  q 


[2M  +  Bnffn4p)  =  5¥  ¥  1  dq  £  Qt (t )  f„ „ (q)  , 


0 


p  i  v  '  ni 


which  is  just  the  Schrodinger  equation  (II-14)  for  the 

2 

Yukawa  potential  with  A  =  G  /4tt 


III. 4  Reduction  to  the  Klein-Gordon  Equation 


Consider  equation  (III-22)  in  the  limit  of 
m2  — >  oo  .  Then, 


2  2 
1  +  P  +  q 

2pq 


=  t 


1 

/  2  ,  2x2 

(q  +  mQ) 


1 

m. 


•  : 


30. 


and 


-  2dni(p2+m|)2  -  -  -  p 


-  [Bni(2mrBnje)  +  p2]  . 


Thus  equation  (ill -22)  becomes , 


pf  +  Bni^2ml  Bn s) 


[P 


G2  1  r°° j  q 


2m 


1 


2m 


1 


]Wp)  “  W  7  •(0dq  p  Vfc)  fni(q)  ’ 


which  is  just  the  Klein-Gordon  equation,,  (11-22),  with 
A  =  G2/4t r  . 


III. 5  Other  Relativistic  Equations 

The  Bethe-Salpeter  equation  for  bound  states 
of  two  identical  scalar  particles  with  mass  m  mediated 
by  a  finite  mass  meson,  is  given  by 


[(i-K+k)2  -  m2][(^K-k)2  -  m2]cpK(k) 


/ aV 

(2  ir)q 


qv(kT)  5 

-§ - ^  '  (III-24; 


(k-k' )  +  [i 


where  p.  is  the  meson  mass. 


G  is  the  coupling  constant. 
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K  is  the  four-momentum  of  the  bound  state  in  the 
centre -of -mass  system  and  k  is  the  relative  momentum. 


5 

Son  and  Sucher  introduced  another  relativistic 
equation  for  a  bound  particle  of  mass  ; 


[2E(p)  -  m.  ]q> (p)  =  -  /  v  (p,  p '  )<p  (p ' )  d3p ' 


(HI-25) 


where 


y(p*p  ) 


n2  n 

G _  m  _ 1 _ m 

(27 r)3  E(p)  C(p-P')2  +  d2]E(p') 


1 

2  2  2 

and  E(p)  =  (p  +  m  )  .  Both  equations  (III-24)  and 

(III-25)  have  the  same  non-relativistic  limit,  the 
equal -mass  Schrodinger  equation  for  a  Yukawa  potential 

v(r)  : 


v(r) 


n2  -ur 
G  e  ^ 

44r  r 


We  shall  use  the  numerical  solutions  of  equations  (ill -24) 
and  (III-25)  given  by  Schwartz^  and  Son-Sucher  respectively, 
to  compare  with  our  results. 


■ 
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Chapter  IV.  METHOD  OF  SOLUTION 


Numerical  solutions  for  the  eigenvalues 
of  equations  (II-14)  and  (II-15)  are  obtained  using 

Q 

Bateman's  method. 


IV. 1  Bateman's  Method 


Consider  the  homogeneous  Fredholm  integral 

equation 


b 

4>(x)  =  A  /  K(x,y)  <t>(y)dy  •  (iV-l) 

a 


To  obtain  the  eigenvalues  A  ,  of  equation  (iv-l), 
Bateman  developed  the  following  procedure.  Define 

b 

T(xrym)  =  f  z)K(z,ym)dz  ,  (lV-2) 

a 


where  i ,m  =  1,2, . n  ,  n  is  a  non-zero  positive 


integer. 

and  x  n  and  y 
i  m 

are  given 

numbers 

in 

the 

interval 

a  <  xrym  <  b  - 

such  that 

T  and 

K 

are 

completely  defined  in  this  interval.  Also  define 


A 


im 


K(x 


+  AT(xrym) 
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Then,  the  eigenvalues  of  equation  (iV-l)  are  given 
approximately  by 

det  (  Aim  ^  =  0  ; 

the  degree  of  approximation  improves  as  n  increases. 

For  the  proof  of  the  above  procedure  and  for  a  discussion 
about  its  range  of  validity,  we  refer  the  reader  to 
reference  (9) • 

IV. 2  Illustration  of  the  Method 


Generally,  in  order  that  the  eigenvalues 
of  an  integral  equation,  such  as  equation  (IV-l),  be 
real  we  must  have  K(x,y)  =  K(y,x)  .  Hence,  we  write 
the  Schrodinger  equation  (II-8)  in  a  symmetric  form 
with,  m,  =  =  m  ,  as 


■2x2 


X^(p)  =  -  (-)  m  /  dp 1 X 0 ( p 1 ) 


v^(p,pT )ppT 


IT 


0 


1  1 
(p2  +  mB) 2 (p  ? 2  +  rnB)2 


(iv-3) 


l 

p  O 

where  X^(p)  =  p(p  +  mB)  u^(p)  .  Note  that  the 

kernel  of  equation  (IV-3)  is  symmetric  since  v^(p,p») 
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is  symmetric,,  and  also  that  equation  (IV-3)  is  of  the 
same  form  as  equation  (iV-l) .  It  turns  out  that  we 
do  not  have  to  symmetrize  our  equation  since  in  this 
method  of  approximation  all  the  unsymmetric  parts  of 
the  kernel  drop  out.  Nevertheless  we  shall  retain 
the  symmetric  form  of  our  equations. 


To  illustrate  the  method  of  approximation 


we  shall  solve  equation  (IV-3)*  with  f  =  0  ,  for 

two  potentials . 

(i)  Delta-function  potential 

Suppose  we  are  given  the  potential 

r* 

v(r)  =  -  A6(r-c)  ,  where  c  is  a  constant.  Then, 

using  equation  (II-10),  for  i=0  ,  we  get 


1 


v0(p,pr)  =  -(|)2^ 


sin(pc)sin(ptc) 

PP' 


Thus,  the  kernel  of  equation  (IV-3)  becomes 


Applying  Bateman Ts  method  to  determine 


the  eigenvalues  A  we  get 
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A= 


20) 


(1  - 


-2o)c  \ 
e  )m 


which  is  the  exact  solution  for  this  potential,  as 

one  can  easily  verify  by  solving  the  Schrodinger  equation 

directly.  From  this  example  we  can  conclude  that 

any  potential  v(r)  for  which  v^(p,pr)  can  be  written 

in  the  separable  form  v^(p,p')  =  A^(p)B^(p')  ,  in  this 

method  of  approximation,  will  yield  exact  solution. 


(ii)  Exponential  potential 

Consider  the  exponential  potential 
v(r)  =  -  A|a2e_^r  .  Equation  (II-10)  then  yields 


v0(p,P’ )  = 


(f>Vs 


2p 


3 


[p2+(p-p' )2] [d2+(p+pr )2] 


and  the  kernel  of  equation  (IV-3)  is 


K(x,y)  = 


4mq~ 


7T 


[i •L2+(x-y)"][^+(x+y)':] 


xy 

2  -i  r  2 


2 


( x2+mB) 2 (y^+mB) 


1 

2 


To  find  the  eigenvalues  A  ,  using  Bateman's  method, 
we  must  solve  det(A^m)  =0  or 


det  { 


[l^+Cx^-yJ2]  [p2+  (x^+ym) 2] 


'i  nr 


3 


00 


W  r  d  _ 

17  0  (q2+mB)  [p2+(x,-q)2]  [p2+(x.+q)2]  [p2+(q-y  )2]  [q2+(q+y  ) 


'i 


. 

■  ? 


3  6. 


o 


(IV-4) 


We  solved  equation  (IV-4)  for  A  with  m  =  p  =  1  , 

^  —  5  and  —  y-^  —  0  ,  X£  —  ho  —  1  ,  .  3  x ^  —  y^  =  4 

The  results  are  tabulated  in  Table  I,  together  with  the 

2 

exact  solution 

=  0 

where  J  (y)  is  the  Bessel  function. 


1°^  bound  state 

2nc^  bound  state 

3rd  bound 

state 

B 

Ae 

Aa 

Ae 

Aa 

Ae 

Aa 

0.00 

1.4458 

1.468 

7.6179 

7.789 

18.7216 

19.228 

0.25 

3.6557 

3.722 

12.3061 

13.045 

1.00 

6.5946 

6 . 658 

17.7115 

18.327 

2.25 

10.1761 

10.277 

4.00 

14.3941 

14.546 

Table  1 .  Value  of  coupling  constant  A  needed  to 
produce  first, second  and  third  S-wave 
bound  state  for  given  binding  energy  B 
(m=[i=l  ;  Ae  -  exact  coupling  constant 
Aa  -  approximate  coupling  constant) 
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As  one  can  see  from  Table  I,  Bateman’s  method  yields 
results  which  are  within  3^  of  the  exact  results.  This 
is  quite  a  good  agreement  considering  the  small  n  (=5) 
we  chose . 

IV. 3  Solutions  of  the  Schrodinger,  Klein-Gordon 
and  Our  Relativistic  Equations 

As  we  have  seen,  equations  (ll-l4)  and  (11-15) 
are  the  Schrodinger  and  the  Klein-Gordon  equation  in 
momentum  space  with  the  Yukawa  potential  v(r)  =  Ae  /r  . 

Now,  in  equation  (II-14),  let  =  mg  =  m  so  that 
2-M  — »  m  .  Also,  we  shall  assume  that  m  in  equation 
(11-15)  represents  the  reduced  mass  of  two  interacting 
particles.  Thus,  when  these  two  particles  have  equal 
mass,  say  m  ,  then  In  equation  (II-15)  2m  — >  m 
Hence,  the  symmetric  form  of  both  equations  (II-14)  and 


Y  f-n) 


2  2 
1+p  +q 

2pq 


00 


xfp)  =  *  ¥  /0d«  V 


(lV-5) 

where  cu  =  mB  for  the  Schrodinger  equation  and  u>  =  B(m-B) 
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for  the  Klein -Gordon  equation.  In  applying  Bateman ;s 
method  to  equation  (IV-5)  we  note  that,  for  Z  =  0  , 
the  kernel  is 


1 


K(x,y)  =  (x2  +  “)2(y2  +  “)2  In1  +  (x+y^ 


1  +  (x-y) 


Thus,  det(A^m)  =  0  implies  that  the  eigenvalues  A 
can  be  obtained  by  solving  the  equation 


,  ,  r  ,  1  +  (xi+yi2 

det  1  In  - 

1  + 


00 


-  / 


dq 


2tt 


In 


1  +  (x^+q) 


0  q2  +  CD  1  +  (x^-q)2 


In 


i  +  (q+ym) 
1  +  (q-ym) 


)  =  0 


(iv-6) 


Similarly,  the  symmetric  form  of  our  relativistic 
equation  (ill -22),  with  mp  =  m2  =  m  and  ^ 
zn^  (p)  =  p(p2+m2)"1//4[d  ^  -  2d„  Jp2+m2)1//2]2  f„,(p)  , 


nZ 


nZ 


nZ 


is 


2  oo 


_  ,  x  ,  2m  r  .  „  /  w  2,  2\  -l/4  /  2  2x-l/4 

zn^  (p)  =  *  ~n~j  Q  dq  Zn,Tq'  +m  )  (q  +m  ) 


r  ,2  n,  /  2  2x  l/2-i  2  r  .,2  0,  ,  2,  2xl/2n2 

*  ^dni  -  2d^(P  +m  )  J  IX*  “  2dM,(q  +m  )  '  ] 


ni 


nZ 


nZ 


n  0  2  ,  0/  2,  2x1/2,  2  2x1/2 

q  [  1  -  2m  +  2(p  +m  )  /  (q  +m  )  j 


2pq 


. 
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where  A  =  G  /bir  .  Hence for  &  =  0  s 


K(x,y)  =  £  (x2+m2)  1'/i,(y2+m2)  -^[d2  -  2d(x2+m2) l/2]  2 


.  [d2  -  2d(y2+m2)1/2]'5  in1  +  (*+?)*  -  [ 

1  +  (x-y)2  -  [  (x2+m2) 1/2  -  (y2+m2)1/2r 


and  det (A^m)  =  0  implies  that 

i  ,  f  ,  n2  r f  2,  2^  1/2  /  2  l/2n2 

1  +  (Xf+3V  "  t(xD+m  )  -  (ym+m  )  ] 


det  {  In 


2  oo 


'i 


,  ,  /  n2  r,  2,  2\  l/2  ,  2,  2sl/2n2 

1  +  "  >-(X0+m  )  "  Urn*111  )  J 


‘4 


mH 


-  A  ^  /  dq(q2+m2)  1//2[d2  -  2d (q2+m2 ) 1//2]  1 


0 


In 


In 


1  +  (x^+q) 2  -  [(x2+m2)1//2  -  (q2+m2) X^2] 2 


IT 


bi 


-  (qW)1'2]2 


1  +  (x^-q)  -  [ (x^+m  ) 

1  +  (q+ymr  -  [(q2+m2)1//2  -  (y^+m2)1^2]2 
1  +  (q-ym)2  -  [(q2+m2)1//2  -  (y2+m2)1//2]2 


}  =  0 


(iv-7) 


where  dQ0  =  d  and  BQq  =  B 

To  find  the  coupling  constant  A  from 
equations  (IV-6)  and  (lV-7),  we  solved  these  equations 
numerically  with  n  =  10  and  x1  =  y1  =  0.01  ,  x2  =  y2  =  1.01  , 

. 3  x^q  =  y^Q  =9*01  .  The  results  are  shown  and 

discussed  in  the  next  section. 
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IV . 4  Results  and  Discussion 

We  chose  our  system  of  units  to  be 

2 

fi  =  c  =  |i  =  1  ,  where  c  is  the  speed  of  light  and 
l_i  is  the  mass  of  the  neutral  meson  (135-01  Mev.)  . 

In  this  system  of  units,  the  mass  of  the  neutron  m 
(939*505  Mev.)  is  6.954  pion  mass,  and  the  binding 
energy  B  is  also  in  pion  mass.  Note  that  we  are  treating 
the  neutron  and  the  meson  as  scalar  particles. 

The  results  of  solving  equations  (IV-6) 
and  (IV-7)  tor  A  ,  with  m  =  6.954  pion  mass,  are 
shown  in  Table  2  and  Table  3  and  are  plotted  in  Fig.  1 
and  Fig.  2  .  These  results  show  that  for  weak  binding 
the  three  equations  agree  with  each  other  well,  whereas 
for  strong  binding  the  difference  is  quite  large.  Note 
that  both  the  Klein-Gordon  and  our  relativistic  equation 
predict  that  for  a  given  A  we  can  have  two  possible 
binding  energies #  This  can  be  explained  by  noting  that 
in  equation  (IV-5)j  if  we  replace  B  by  m-B  the 
equation  remains  unchanged.  Thus  A  is  symmetric 
about  B  =  m/2  for  the  Klein-Gordon  equation,  and  for 
a  given  A  both  B  and  m-B  are  acceptable  binding 
energies.  Now,  since  our  relativistic  equation  reduces 
to  the  Klein-Gordon  equation  in  the  limit  that  mg  is 


. 


' 

- 


large,  we  expect  it  to  retain,  at  least  partially, 
this  peculiar  behaviour  of  the  Klein-C-ordon  equation 
even  for  rrig  = 

Table  4  shows  the  variation  of  A  with 
l/m  ,  in  the  ground  state,  for  different  binding 
energies  B  .  In  Fig.  3  we  plotted  these  results. 
From  Fig.  3  we  note  that  as  m  tends  to  infinity, 
the  three  equations  agree  with  each  other  very  well 
for  all  B  .  For  m  finite  (but  not  less  than  one 
pi on  mass),  however,  the  agreement  is  good  for  weak 
binding  (B  smaller  than  one  pion  mass)  and  very  poor 
for  strong  binding  (B  larger  than  one  pion  mass) . 

Finally,  to  compare  our  results  with  those 
obtained  by  Schwartz  and  by  Son  and  Sucher  we  solved 
our  equations  with  m  =  *  l  .  The  results  are 

tabulated  in  Table  5  and  plotted  in  Fig.  4  . 
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B 

AS 

ak-g 

ar 

0.00 

0.249 

0.249 

0.248 

0.01 

0.326 

0.326 

0.323 

0.04 

0.406 

0.406 

0.402 

0.0  9 

0.489 

0.487 

0.483 

0.16 

0.569 

0.565 

0.559 

0.25 

0.649 

0.642 

O.616 

0.36 

0.726 

0.713 

0.700 

0.49 

0.803 

0.785 

0.771 

0.64 

0.884 

0.855 

0.837 

0.81 

0.964 

0.920 

0.901 

1.00 

1.042 

0.983 

O.961 

1.21 

1.120 

1.042 

1.020 

1.44 

1.205 

1.089 

I.O78 

1.69 

1.287 

1.148 

1.124 

.  1.96 

1.372 

1.194 

1.179 

2.25 

1.455 

1.234 

1.224 

2.56 

1.541 

1.264 

1.262 

2.89 

1.631 

1.287 

1.297 

3.24 

1.728 

1.299 

1.331 

B 

AS 

. -  ■  ■  - 

ak-g. 

A  R 

3.61 

1.801 

1.301 

1.354 

4.00 

1.896 

1.290 

1 . 366 

4 .41 

1.989 

1.262 

1.376 

4.84 

2.086 

1.216 

1.379 

5.29 

2.184 

1.143 

1.375 

5.76 

2.284 

1.038 

1.361 

6.25 

2.385 

0.880 

1.335 

6.76 

2.487 

0.595 

1.297 

7.29. 

2.590 

1.258 

7.84 

2.697 

1.208 

8.4l 

2.791 

1.174 

9.00 

2.901 

1.062 

9.61 

3.029 

0.968 

10.24 

3.144 

0.832 

10.89 

3.260 

0.732 

11.56 

3.378 

0.592 

12.25 

3.498 

o.46o 

12.96 

3.622 

0.269 

13.69 

3.747 

Table  2.  Coupling  constant  A  required  to  produce 
the  first  S-wave  bound  state  of  two  neutrons  of  mass  m 
(=  6.954  pion  mass)  with  given  binding  energy  B  in 
pion  mass.  (Ag  -  from  Schrodinger  equation;  AR_G  -  from  the 
Klein-Gordon  equation;  AR  -  from  our  relativistic  equation) 


' 


43. 


Fig.  1.  Ground  state  coupling  constant  A 
versus  the  binding  energy  B  in  pion  mass  unit. 

SE  -  from  the  Schrodinger  equation 
RE  -  from  our  relativistic  equation 
K-G  -  from  the  Klein-Gordon  equation 
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B 

As 

ak-g 

"ar 

0.00 

0.944 

0.944 

0.955 

0.01 

1.132 

1.132 

1.135 

0.04 

1.336 

1.334 

1.359 

0.09 

1.485 

1.482 

1.514 

0.16 

1.643 

1.634 

1.678 

o 

• 

ro 

Ul 

1.852 

1.835 

1.804 

0.36 

2.041 

2.010 

2.052 

0.49 

2.239 

2.189 

2.261 

0.81 

2.651 

2.536 

2.640 

1.00 

2.896 

2.703 

2.874 

1.21 

3.131 

2.895 

3.096 

1.44 

3.346 

3.101 

3.283 

1.96 

3.791 

3.315 

3.536 

2.25 

4.094 

3.425 

3.716 

3-24 

4.898 

3.537 

4.242 

4.00 

5.366 

3.507 

4.456 

4 .4l 

5.749 

3.438 

4.542 

4.84 

6.043 

3.377 

4.597 

5.29 

6.388 

3.179 

4.623 

Table  3.  Coupling  constant  A  required  to 

produce  the  second  S-wave  bound  state  of  two  neutrons  of 

mass  m  (=  6. 954  pion  mass)  with  given  binding  energy  B 

in  pion  mass.  (Ac  -  from  the  Schrodinger  equation;  n  - 

from  the  Klein-Gordon  equation;  A0  -  from  our  relativistic 

n 

equation . ) 


Fig.  2.  First  excited  state  coupling 
constant  A  versus  the  binding  energy  B  in  pion 
mass  unit. 

SE  -  from  the  Schrodinger  equation 
RE  -  from  our  relativistic  equation 
K-G  -  from  the  Klein-Gordon  equation 
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Fig.  3. 

of  the  mass,  l/m  ,  in  (6.95^  pion  mass)-'1'  for  given 
binding  energy  B  in  pion  mass  unit. 

SE  -  from  the  Schrodinger  equation 
K-G  -  from  the  Klein-Gordon  equation 
RE  -  from  our  relativistic  equation 
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Fig.  4.  Ground  state  coupling  constant  A 
versus  the  binding  energy  B  for  m  =  p.  =  1  . 

SE  -  from  the  Schrodinger  equation 
RE  -  from  our  relativistic  equation 
K-G  -  from  the  Klein-Gordon  equation 
B-S  -  from  the  Bethe-Salpeter  equation 
S-S  -  from  the  Son-Sucher  equation 
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Chapter  V.  CONCLUSIONS 


As  we  have  already  seen,  the  values  of  the 
coupling  constants,  for  weak  binding  (B  «  m)  ,  predicted 
by  the.  lowest -order  "N-quantum  approximation"  are  in 
good  agreement  with  those  predicted  by  the  Schrodinger 
equation.  On  the  other  hand,  as  was  pointed  out  by 
Vosko^,  the  Bethe-Salpeter  equation,  for  very  small 
values  of  B  , yields  values  for  the  coupling  constants 
which  are  about  40^?  larger  than  those  obtained  from 
the  Schrodinger  equation. 

For  strong  binding  (B  >  m) ,  the  agreement 
of  all  the  relativistic  equations  with  the  Schrodinger 
equation  is  poor.  This  is  not  surprising  since,  a's 
was  pointed  out  by  Greeberg^ ,  a  particle  which  is. 
strongly  bound  by  a  Yukawa  potential  must  move 
relativistically  and  thus,  must  be  described  by  a 
relativistic  equation.  Moreover,  since  all  of  the 
relativistic  equations  reduce  to  a  Yukawa-like  interaction 
in  the  proper  limit,  the  particles  described  by  these 
equations  must  move  relativistically  and  independently 


of  m  . 


* 
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To  see  if  one  can  represent  the  interaction 
predicted  by  the  "N-quantum  approximation" ,  in  terms  of 
a  potential,  we  write  both  the  Schrodinger  and  the 
Klein-Gordon  equations  in  the  most  general  (non-local) 
form  as  (see  Appendix) 

2  _ 

(  m  +  03  )u(p)  =  '  /d^q  v(P*q)  u(q) 

where  a)  =  b  for  the  Schrodinger  equation  and 

=  B(m-B)/m  for  the  Klein-Gordon  equation.  We  also 
rewrite  our  relativistic  equation  (III-21)  as 

(£2+  m) f (7)  -  .  2GV  /d3q  -  (P2Au.iifl3i - -  , 

(2tt)  d  qQ(d  -2dp0)  (l-2m  +2pQq0-2p.q) 

where  pQ  =  (p  +m  )  '  and  qQ  =  (q  +m  )  ' 

Comparing  the  above  two  equations  we  get 

v (r,  r 1 )  =  (2tt)  ~3//d3pd3q  e~ip*rv(p,q)eiq'r 

cd)  eiq 

+2pQq0-2p.q) 

To  approximate  the  above  integral  we  expanded  the  kernel, 

o 

v (p, q) j  in  terms  of  m  (see  Appendix).  Thus,  up  to  l/mJ 


(2ir)  3//d3pd3q 


e-ip-r(p2/m  + 

qQ(d2-2dpn) (l-2m2 


'0- 


- 
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we  have 


Jv(r,rx  )¥(r*  )d3rr  =  -  ^  (l  +  ^  +  — 2  +  — . .)  —  ¥(r)  +  F-.'F(r) 

/  !  m  ,W  wi  — ' 


B  .  B 


B3  ,e'r 


4m  8m~ 


for  the  Schrodinger  equation,  and 


/v(r,r 1 )  f (r 1 ) d3r 1  =. - 


2  2 
G  (1  _  B  B 


B3  ,e'r 


4ur 


2m  4m2 


8m 


■3)V  +  F2f^ 


for  the  Klein-Gordon  equation,,  where  (see  Appendix)  F-^ 
and  Fg  are  complicated  functions  of  B,  m/  r;  and 
p  (=  -iV) .  Thus,  we  see  that  the  lowest-order  N-quantum 
approximation  predicts  a  non-local  energy-dependent 
interaction.  For  m  very  large  the  non-locality  tends 
to  zero  and  up  to  order  of  l/m  we  have 


v(r,rT ) 


(1  + 


-r 

—  5 (r-r 1 ) 
r  v  ' 


for  the  Schrodinger  equation  and 


G2 

4¥ 


v(r,r' )  =  -  ^  (1  - 


-r 


6 (r-r f ) 


for  the  Klein-Gordon  equation.  Hence,  l/m  roughly 
represents  the  range  of  the  non-locality  of  the 


interaction . 


;  '  ii 


Finally,  we  note  that  one  can  use  equation 

(III-21)  and  Bateman’s  method  to  determine  numerically 

the  wave  functions  of  the  bound  states.  This  in  turn 

would  help  us  to  understand,  at  least  qualitatively, 

the  overall  nature  of  the  interaction  (repulsive  or 

attractive)  and  provide  an  interesting  comparison 

to  the  wave  functions  one  obtains  from  the  Bethe- 

12,13 


Salpeter  equation 
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APPENDIX 


THE  FORM  OF  THE  INTERACTION  PREDICTED 
BY  OUR  RELATIVISTIC  EQUATION 


The  Schrodinger  and  the  Klein-Gordon  equations 
with  the  non-local  potential  v(r,r')  can  be  written 
as 

2 

(-  ^  +  <^>)L(r)  =  -  /  d^r’  v (r, I?1 )  f  (r 1 )  (A-l) 

where  =  B  for  the  Schrodinger  equation  and 
=  B(m-B)/m  for  the  Klein-Gordon  equation.  Let 

v(r,rt)  =  (2tt)  "^//d^pd^q  e  ip‘r  vf^q)  e1^'1  (A-2) 

=  T[v(p,q) ]  . 

Then, 

v(p4q)  =  (2tt)  _3//d3rd3rl  eip”r  v(r,vx)  e  (A-3) 

=  T’^vp^r' )  ]  . 
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Also  let 

¥(r)  =  (2t r)~3//2  f  d3p  ufp)  e"ip,r  .  (A-4) 


Substituting  equations  (A-4)  and  (A-2)  into  equation 
(A-l)  we  get 


(  £  +  cw  )u(p)  =  -  /  d3q  v(p,q)  ufq)  .  (A-5) 


Consider  now  equation  (ill -21)  with  rrr  =  m 2 


(d2  -  2dp0)f(p)  = 


where  d  =  2m  -  B 

and  fni®  =  fS) 
as 


2Gjn2  /d3q _ IS) _ 

(2tt)3  qQ(l  -  2m2  +  2pQq0  -  2p.q) 

P0  =  (p2  +  m^)1^2  ,  qQ  =  (q2  +  m2) 
This  equation  can  be  rewritten 


(f  +  *)?(?)  ,  -  2Gfm2  /d3q  - (}>?/*+.  - 

(27 r)J  qQ(d  -2dp0) (l-2m  +2pQq0-2p.q) 

(A -6) 

Comparing  equation  (A -6)  with  equation  (A-5)  we  see 
that 


v(pS) 


G2  m^_  _ (p2/m  +  a)) _ 

^  tt2  q0(d2-2dpQ)  (l-2m2+2p0q0-2p.q) 


1/2 


(A-7) 


56. 


and  using  equation  (A-2)  we  have 


v(r3r' ) 


2  2 

G  m  So^-3 

7 r 


(2tt)  "^//d^pd^q 


e-1?~?  fp2/m  +  oQ  e1^7' 

q0(d2-2dp0) (l-2m2+2p0q0-2p.q) 


(A-8) 


We  proceed  to  approximate  equation  (A-8)  by  assuming 
that  m  is  large  and  thus  expanding  the  kernel,  v(p,q), 
in  powers  of  l/m  .  Hence, 


w  /  2  ,  2x1/2 

l/q0  =  (q  +  m  ) 


-  [ 
m 


1  + 


a-,  ap 

-2  +  -zr 

m  m 


) 


2  4 

where  a1  =  -  q  /2  and  a2  =  3q  /8  ; 

[1  -2m2+2  (p2+m2) 1/'2  (q2+m2)  1^2-2p .  q]  _1 


=  [1  +  (?-cf)2 


-1 


[1  + 


bl 


+ 


where 


,  2  2  s  2 

b,  =  — ^ ‘ - p—  and  b 

1  4 [ 1  +  (?-q)2] 


(P2-^4  tpW- 

8[l+(jr-q)  ] 


,2  2,2 
(p  -q  )  } 

2[l+(p-5)2] 


and 


(d2  -  2dp0)  1  =  (2m-B)  1[2m-B-2(p2+ffi2)1//2]  1 


■ 
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-  (2mB)  [  1  +  +  —?2  i-  — r  f  — 7T 

m  rrr  m 


ci  .  c2  ^3  c4  ] 

+  Q  "r  JT  -I  } 


where  c1  =  -  p2/B  +  B/2  ,  c2  =  B2/4  +  p^/B2  -  p2/2  , 
3  ;q  ,  J^ot,  ,  ^2,)t  6  .3 


=  BJ/8  +  p  /2B  +  p  /4B  -  Bp  /4  -  p°/Bj  and 
c4  =  B4/l6  -  p°/2B2  +  p8/B4  +  p4/8  -  p^/2B2  +  p4/4  -  B2p2/8 
Hence,  up  to  l/mJ  v(p,q)  can  be  written 


as 


G2  1  1  1 

v(p,q)  =  -  WZ 2  2B 

7T 


V1  V?  VQ 

^2  tv0  +  TT  +  T  +  3]  -• 

±+(p-qj  m  rrr 


2 


2  2  2 

where  v„  =  B  ,  v,  =  B2/2  ,  vQ  =  B3/4  +  tt[  ^P  ~q  4  -  2qc] 

1  2  4  l+(p-q)2 


4  4  2  ,22\2 

and  Vo  =  §  +5  +  3  [  2  ~  2c^  for 

J  ‘  1+ ( P-q) 

Schrodinger  equation  and 
r0  -  -  *  vi 


3 


vn  =  B  ,  v,  =  ~  B2/2  ,  v,  =  -  +  |[  (p2~q^  -  2q2  +  4p2]  , 

l+(p-q) 


and  -  -  b4  3p4  fi2[  (p2-32 

^  V3  -  S  T  -5'  1+(M)2 


-  2q  ]  for  the 


Klein-Gordon  equation.  Now,  using  the  identities 


T  [- 


-r 


2t r  l+(p-q) 


■2 


]  =  tt-  6  (r-r 1 ) 


and 


T  [ 


0  - - — p-p  ]  =  e  r5  (r-r 1 ) 

tt2  [  1+  (p"-q) 2  ] 2 
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we  get 


v(r,r* )  =  - 


2  2  R  —  p 

G  f  n  B  ,  B  n  e  x,  \ 

^  1  (1  +  2-  +  ^  +  ^3}-  6^-r  > 


+  ^  V^,[5(r-r')l  -  ^3  (?-?')  ] 


+  -dy(l  +  fm)(VB  -  V2, ) 2 [e'r6 (?-?')]  3 


4m 


for  the  Schrodinger  equation  and 


v(r,r  )  =  - 


2  2 
G  {  (l  -  B  B 


Wir 


2m  2 


4in  8m 


R3  e~r  -> 


3/  r 


+  -  lm)(Vr  "  (*-?■)]  +  Ay(l  - 


8m 


2m“ 


2m- 


+  (— ■ ,  d  -  f,  V^)[t-r6(r-r')l  ) 


4Bm 


3  r  2  v'  u  r 


m 


8 


for  the  Klein-Gordon  equation  with  V  =  and  V^r 


Thus,  in  equation  (A-l) 


fd^r'v(r,r'  )f  (r! )  -  -  ^(1  +  |m  +  ~2  +  + 


4m  8m~ 


for  the  Schrodinger  equation,  and 


•  v;t [5(r-r! ) ] 


3 

c^F1  * 


Fni(r)  , 
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/ d^r J  v(r, r * ) IF (r ' ) 


2  2  ?  -r 

G  / -.  B  B  B\6  .t,  /— ►n  t-,  ,T/  / 

-  -  2m  -  ^2  -  ^)~  +  VW 


for  the  Kiein-Gordon  equation  where 


1  .  B  \  e  r‘  2  B  „4  r  e  r 


F,  =  -^g(l+f  )%-  V£  -  -^W[S-  ]  +  -~(1+|  )  (V4e  r-2V2e'rV2+e‘rV4) 
2m2  2ra  r  4nP  r  8rrt  2m 


and 


1  /■.  B  se  r„2  3  „4re  rn  ,  1  B  w„4„-r  on2„-r„2,,-rJ1 


F2  =  -^-aJ-r  v 

2  2m  4m- 


(V  e  -2V  e  V  +e 

J  r  8m2  2m 


and  thus  we  see  that  the  interaction  is  non-local  and 


energy-dependent . 


.. 


